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摘要 ”为 了 研究 过 渡 流 下 后 向 台阶 的 流动 与 瞬 态 传 热 特性 ， 建 立 了 二 维 不 可 压缩 后 向 台阶 流动 的 数值 传 热 模型 ， 采 用 非 
定常 数值 方式 模拟 了 后 向 台阶 的 内 部 流 场 , 比较 了 台阶 下 游 底 面 瞬时 努 塞 尔 数 Nu 和 时 均 努 塞 尔 数 Nu 的 分 布 , 分 析 了 台 
阶 下 游 通道 内 温度 分 布 规律 。 结 果 表 明 : 过 渡 流 下 台阶 底面 Nu 数 沿 着 流动 方向 出 现 多 个 峰值 ，Nu 数 始终 围绕 Nu 数 发 
生 波动 。Re 数 为 700 时 ， 主 回流 区 未 出 现 二 次 回流 ，Nu 数 在 主 回流 区 下 游 才 发 生 波动 ，Re 数 为 1200 时 ， 主 回流 区 出 现 
二 次 回流 ，Nu 数 在 主 回流 区 就 发 生 了 波动 。 台 阶 下 游 通道 内 温度 沿 法 向 呈 阶 梯 式 变化 ，Re 数 越 大 ， 近 底面 处 的 温度 变化 
梯度 也 越 大 。 
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Study on Flow Field and lnstantaneous Heat Transfer Characteristics for 


Backward-facing Step in Transition Flow 

ZHONG Min-Bo” XIE Wei-An YUANYin-Nan'” XIGuan-Nan'” 
(1. School of Energy and Power Engineering, Jiangsu University, Zhenjiang 212013, China; 
2. School of Mechanical Engineering, Nantong University, Nantong 226019, China) 
Abstract In order to study the flow field and the instantaneous heat transfer characteristics for 
backward-facing step in transition flow, a two-dimensional incompressible numerical model of 
backward-facing step flow was established. The internal flow field was studied with the unsteady numerical 
method. The instantaneous Nusselt number (Nu) and the time-averaged Nusselt number (Nunm) distributed on 
the bottom wall were compared. The temperature distribution of the channel was analyzed. The results show 
that Nu appears multiple peaks on the bottom wall along the flow direction in transition flow. The distribution 
of Nu always fluctuates around the Num. When the Reynolds number is 700, Nu fluctuates only in the 
downstream of the primary recirculation region because there is no secondary recirculation zone in the primary 
recirculation region. When the Reynolds number is 1200, Nu fluctuation occurs in the primary recirculation 
region because secondary reflux appears in the primary recirculation region. The temperature changes in the 
form of step in the normal direction for the channel downstream of the step. The larger the Reynolds number is, 
the greater the temperature gradient 1s. 
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0 引言 


后 向 台阶 流动 分 离 现 象 普遍 存在 于 动力 、 化 工 、 数 Re 增 大 ， 流 场 具 有 三 维特 性 。 之 后 ， 众 多 科研 工 
石油 与 能 源 等 领域 。 迄 今 ， 针 对 后 向 台阶 的 流动 特 ”作者 对 后 向 台阶 的 流动 特性 进行 了 更 加 深入 细致 的 
性 已 有 大 量 的 实验 和 数值 研究 ，Armaly 等 由 通过 实 。 ”研究 。Williams 和 BakerD 对 层 流 状态 下 的 流 场 进 行 


验 和 数值 方式 首次 对 同一 模型 研究 后 发 现 随 着 雷诺 了 三 维 数值 模拟 ，Tylli 等 中 研究 了 台阶 侧 壁 对 流 场 
收 稿 日 期 : 

: 国家 自然 科 和 起 金 资助 项 目 (No.51176079) 

作者 简介 : 仲 敏 波 (1981-)， 女 ， 讲 师 ， 博 士 生 ， 主 要 从 事 传 热 强化 研究 。 通 讯 作 者 : 喜 冠 南 ， 教 授 ，guannanxi@ntu.edu.cn 


4 影响 ; Erturk 四 通过 数值 模拟 研究 了 Re 数 变化 对 流 
体 再 附着 点 位 置 的 影响 。 国 内 学 者 王 兵 等 牛 模 拟 了 满 
流 时 后 向 台阶 流 场 中 大 尺度 涡 结构 的 发 展 和 演化 过 
程 ， 齐 鄂 荣 等 "应 用 PIV 技术 对 层 流 、 过 渡 流 、 注 
流 状 态 下 流 场 中 旋涡 的 生成 、 发 展 和 演化 规律 进行 
了 研究 。 

但 是 到 目前 为 止 ， 只 有 少数 研究 是 针对 后 向 台 
阶 的 传 热 特性 ， 且 研究 范围 主要 限于 层 流 I 和 灌流 
0 状态 。 然 而， 过 渡 流 作 为 一 种 特殊 的 流动 状态 
其 流动 和 传 热 特性 有 别 于 层 流 和 满 流 4， 所 以 研 
究 过 渡 流 状态 下 的 流动 传 热 问题 对 理论 研究 及 工程 
应 用 都 具有 重要 的 意义 。 

本 文 根 据 过 渡 流 的 非 定 


党 特点， 建立 了 二 维 非 
定常 流传 热 数值 计算 模型 ， 模 拟 了 过 渡 流 下 的 流 场 
和 温度 场 , 对 台阶 底面 暗 时 努 塞 尔 数 UNu) 和 时 均 努 大 
尔 数 (Num) 的 分 布 作 了 对 比 ， 分 析 了 台阶 下 游 通道 内 
的 温度 分 布 规律 。 

1 计算 模型 和 计算 方法 

1.1 计算 模型 和 边界 条 件 

图 1 表示 后 向 台阶 流动 的 数值 计算 模型 。 坐 标 
原点 位 于 台阶 垂直 壁面 和 下 游 壁面 的 交点 处 ， 治 着 
流体 流动 方向 为 x 方向 ， 垂 直 于 流体 流动 方向 为 y 
方向 。 台 阶 高 度 为 S， 平 板 间 高 度 为 互 ， 台 阶 扩张 比 
ER= 互 /( 五 -=S)=2 。 人 台阶 上 游 进口 长 度 厂 =0.5H ， 
台阶 下 游 出 口 长 度 L =50H 。 


图 1 计算 模型 
Fig.1 Calculation region 

台阶 进口 边界 温度 均匀 分 布 T=10'C, 进口 速度 
uin 充分 发 展 ， 满 足 抛物 线 分 布 ， 出 口 边界 的 速度 场 
和 温度 场 符合 边界 层 近 似 理论 中， 壁面 边界 : 台阶 
所 有 壁面 速度 无 滑 移 ， 台 阶 下 游 底面 温度 恒定 
Tw=40C， 其 余 壁 面 满 足 绝热 边界 条 件 。 
1. 2 控制 方程 及 求解 


以 常 物性 二 维 连续 不 可 压缩 非 稳 态 流体 (空气 ) 
为 研究 对 象 ， 采 用 的 控制 方程 包括 质量 守恒 、 动 量 
守恒 、 能 量 守 恒 方 程 ， 表 达 如 下 : 

6 和 
质量 守恒 二 EN (1) 
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上 述 控制 方程 (1)~(4) 中 ，u、v 分 别 代表 x、y 方 
向 速度 分 量 , P 是 压力 , T 是 温度 ,t 表示 时 间 ，p、 
HU、A、 Gp 分 别 代表 流体 密度 、 烙 性 系数 、 导 热 系 数 、 
定 压 比 热 。 

采用 有 限 容积 法 求解 二 维 非 定常 控制 方程 ， 动 
量 方程 中 对 流 项 和 扩散 项 分 别 采用 三 阶 迎 风 差 分 
(QUICK 算法 ) 和 二 阶 中 心 差分 求解 ， 时 间 项 采用 隐 
式 求解 法 。 求 解 方程 时 , 每 一 个 时 间 步 长 中 采用 ADI 
算法 ， 和 迭代 过 程 中 ， 压 力 和 速度 的 修正 均 采 用 
SIMPLE 算法 。 通 过 课题 组 编制 的 FORTRAN 程序 实 
现 控制 方程 组 的 求解 。 
1.3 网 格 划分 及 模型 验证 

采用 非 均 匀 网 格 划 分 计算 区 域 ， 
化 梯度 大 的 台阶 分 离 点 和 上 、 下 壁面 附近 区 域 对 网 
格 加 密 。 为 了 消除 网 格 差异 对 计算 精度 的 影响 ， 进 
行 了 网 格 独立 性 计算 。 网 格雷 诺 数 Res 反映 了 网 格 的 
疏 密 情况 ， 考 虑 到 计算 机 的 容量 和 计算 速度 ， 经 过 
计算 确定 在 整个 研究 范围 内 Res 取 值 12。 采 用 这 
网 格 精度 计算 得 到 的 再 附着 点 位 置 与 他 人 的 数值 和 
实验 结果 对 比如 图 2 所 示 ， 比 较 结 果 表 示 目 前 的 计 
算 方法 是 正确 可 靠 的 。 


在 流动 参数 变 
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图 2 再 附着 点 位 置 比较 


Fig.2 Comparison of reattachment point positions 
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图 3 
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Fig.3 Vortex field 


2 计算 结果 与 分 析 


本 文 研究 了 Re 数 ( 以 台阶 高 度 $ 作为 特征 长 度 ; 
速度 取 进 口 平 均 速度 元 ) 500、700、1200 三 种 不 同 
工 况 下 后 向 台阶 的 流 场 和 传 热 特性 。 图 3 表示 这 三 
种 工 况 下 的 瞬 态 涡 量 场 。 从 图 上 可 以 看 出 ，Re 数 发 
生变 化 ， 流 场 的 状态 也 发 生 了 变化 。Re=500 时 ， 人 台 
阶 下 游 未 产生 旋涡 ， 说 明 流 场 处 于 层 流 状态 ; 
Re=700、1200 时 , 台阶 下 游 产 生 了 交错 旋转 的 旋涡 ， 
说 明 流 场 处 于 过 渡 流 状态 。 为 了 比较 后 向 台阶 在 过 
渡 流 与 层 流 状态 下 的 瞬 态 传 热 特性 ， 文 中 对 Re=500 
时 后 向 台阶 的 流 场 和 有 瞬 态 传 热 特性 规律 同时 进行 了 
研究 。 图 4 表示 与 图 3 对 应 时 刻下 的 速度 矢量 场 ， 
图 中 的 黑 线 表示 流 线 ,“ 全 ”表示 瞬时 再 附着 点 位 置 。 
从 图 上 可 以 看 到 台阶 下 游 流 场 〈 旋 涡 发 展 ) 随 着 Re 
数 增 大 的 发 展 过 程 。 图 中 Re=500 的 层 流 状态 下 ， 后 
向 台阶 下 游 底面 只 在 靠近 台阶 处 存在 一 个 主 回流 
区 ， 再 附着 点 下 游 流 场 重新 发 展 ，Re 数 增 大 后 ， 
Re=700， 主 回流 区 内 涡 核 发 生 分 裂 ， 同 时 再 附着 点 
下 游 流 场 沿 着 流动 方向 产生 多 个 旋涡 Re 数 进一步 
增 大 ，Re=1200， 主 回流 区 产生 了 二 次 回流 区 ， 同 时 
在 台阶 角落 xw/H 志 1 的 区 域 有 角 涡 形成 。 

图 5 表示 台阶 下 游 底 面 时 均 努 塞 尔 数 Nun 沿 着 
流动 方向 的 分 布 。 从 图 上 可 以 看 出 ，Re=500 时 ，Nunm 
数 沿 着 流动 方向 只 有 一 个 峰值 ， 符 合 层 流 时 的 分 布 
特征 。Re=700、1200 时 ，Nun 数 沿 着 流动 方向 出 现 
两 个 峰值 (上游 的 称 为 主峰 、 下 游 的 称 为 次 峰 )， 主 
峰值 大 于 次 峰值 。Re 数 由 700 增加 到 1200， 主 峰值 
由 5.9 增加 到 10.3， 次 峰值 由 5.0 增加 到 8.3， 主 峰 
值 的 增加 程度 是 次 峰值 增加 程度 的 1.33 倍 。 主 峰 位 
置 由 x/H=8.0 变化 到 xw/H=6.6， 次 峰 位 置 由 WE=13.1 


~ 


变化 到 x/H=9.7， 两 峰 位 置 随 着 Re 数 增 大 都 向 着 台 
阶 方向 移动 ， 但 次 峰 的 移动 速度 是 主峰 移动 速度 的 
2.43 倍 ， 使 两 峰之 间距 离 随 着 Re 数 的 增加 而 减 小 。 


wH 
(c) Re=1200 


图 4 流 线 和 速度 矢量 图 


Fig.4 Streamline and the velocity vector 


Re=500 
~--~- Re=700 


图 5 台阶 底面 时 均 努 塞 尔 数 Nu 分 布 


Fig.S Time-averaged Nun distribution on the bottom wall 
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图 6 表示 台阶 下 游 底 面 时 均 摩擦 系数 Cm 沿 着 流 
动 方向 的 分 布 。 图 中 Cpw<0 的 区 域 是 回流 区 。 通 过 
Cn 的 分 布 可 以 看 到 : Re=500 时 , 台阶 底面 只 存在 一 
个 主 回流 区 ;Re=700 时 ， 主 回流 区 下 游 出 现 了 一 个 
二 次 回流 区 ; Re=1200 时 , 主 回流 区 下 游 的 二 次 回流 
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态 流 场 的 时 刻 。 从 图 中 可 以 看 出 ，Re=500 的 层 流 状 
态 下 ，Nu 数 和 Nun 数 沿 着 流动 方向 的 分 布 重合 。 而 
在 过 渡 流 状态 下 ，Nu 数 沿 着 流动 方向 会 发 生 剧烈 的 
波动 ， 出 现 多 个 峰值 : Re=700 时 ，Nu 数 在 主 回流 区 
先 快速 增 大 ， 在 主 回流 区 下 游 才 发 生 波 动 ; Re=1200 


区 消失 ， 但 主 回流 


区 内 部 出 现 了 二 次 


回流 区 。 所 以 ， 


时 ，Nu 数 在 主 回流 


区 就 发 生 波动 ， 沿 着 流动 方向 波 


从 时 均 的 角度 研究 二 次 回流 区 ， 随 着 Re 数 的 增 大 ， 


动 振 幅 快 速 增 大 ， 在 有 


附着 点 附近 到 达 最 大 ， 然 后 


台阶 底面 的 


次 回流 区 将 由 主 回 


流 区 下 游 向 主 


区 内 部 迁移 。 


回流 


因此 ， 图 5 中 主峰 、 次 峰 位 置 随 着 Re 


数 的 增 大 向 台阶 方向 移动 是 二 次 回流 区 向 上 游 迁移 


又 逐渐 减 小 。 图 4(b) 中 主 回流 段 未 产生 二 次 回流 , 对 
应 图 7(b) 中 该 瞬时 的 Nu 数 在 回流 段 没有 发 生 波动 ， 
与 Nun 的 分 布 重合 .图 4(c) 中 回流 段 底 面 出 现 了 三 个 


后 


的 结果 ， 而 二 次 


回流 区 出 现在 主 回流 区 内 对 主峰 


附 壁 涡 ， 对 应 图 7(c) 中 该 瞬时 的 Nu 数 在 该 


立 置 出 现 


值 的 快速 增加 又 起 到 了 促进 作用 。 
图 7 表示 台阶 下 游 底面 瞬时 努 塞 尔 数 Nu 治 着 流 
动 方向 的 分 布 。 图 中 的 瞬时 Nu 数 分 布 对 应 图 3 中 瞬 


了 三 个 波峰 ;图 4(c) 中 瞬时 再 附着 点 位 于 WH=7.8 处 ， 
图 7(c) 中 该 位 置 的 Nu 数 也 到 达 峰 值 。 计 算 结 果 显示 : 
Nu 数 的 波动 始终 围绕 着 Nun 数 的 分 布 。 


Re=500 
--~ Re=700 
--—--Re=1200 


台阶 底 国 


图 6 


i 时 均 摩擦 系数 Cm 分 布 


图 8 表示 台阶 下 游 底 面 瞬时 摩擦 系数 Cr 沿 着 流 
动 方向 的 分 布 。 从 图 上 可 以 看 到 在 某 一 瞬时 台阶 底 
面 回流 区 的 分 布 情况 : Re=500 时 ， 流 场 的 瞬时 分 布 
和 时 均 分 布 表现 一 致 ， 说 明 层 流 时 流 场 比较 稳定 ; 
Re=700、1200 的 过 湾流 态 下 ， 沿 着 流动 方向 后 向 台 
阶 下 游 底面 存在 多 个 回流 区 ， 这 一 现象 与 流 场 的 时 
均 表 现 不 同 ， 说 明 过 渡 流 下 ， 流 体 流 动 状态 复杂 ， 
流 场 的 影响 因素 多 。 不 同 流 态 下 ，Cf 数 与 Cm 数 分 布 
之 间 的 关系 和 Nu 数 与 Nun 数 分 布 之 间 的 关系 类 似 。 
说 明了 流动 与 传 热 之 间 存 在 密切 关系 。 


Fig.6 Time-averaged Cm distribution on the bottom Wall 
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图 8 


Fig.8 Instantaneous Cr distribution on the bottom wall 
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图 9 温度 分 布 云图 


Fig.9 Contour of temperature 
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图 10 温度 


在 法 向 的 分 布 


yc WO ye 


Fig.10 Distribution of temperature in the normal direction 


为 了 对 后 向 台阶 流动 的 传 热 特性 有 更 加 全 面 的 
认识 。 图 9 给 出 了 后 向 台阶 通道 内 的 温度 分 布 云图 ， 
此 瞬 态 温度 分 布 时 刻 对 应 图 3 中 的 瞬 态 流 场 时 刻 。 
比较 图 3 和 图 9, 由 于 流 场 结构 会 影响 温度 场 的 分 布 ， 
温度 场 的 分 布 特征 和 流 场 的 分 布 特征 类 似 。Re=500 
时 ， 流 场 处 于 层 流 状态 ， 温 度 分 布 呈 明 显 的 层 状 。 


图 上 可 见 ， 通 道 
式 ， 在 不 同 的 区 : 


内 温度 在 法 向 的 变化 呈 明 显 的 阶梯 
成 ， 各 截面 上 温度 在 法 向 的 变化 也 


各 不 相同 。 
发 生变 化 ， 人 台阶 


度 沿 法 向 


而 在 过 渡 流 时 ， 由 于 流 场 中 的 旋涡 运动 引起 台阶 底 
和 的 热流 体 和 上 部 冷 流体 混合 ， 流 道内 温度 分 布 出 
现 波动 ， 在 近 壁 面 处 温度 梯度 最 大 。 在 再 附着 点 上 
游 ， 温 度 分 布 表 现 为 两 种 形式 : Re=700 时 ， 温 度 在 
台阶 高 度 以 下 通道 呈 层 状 分 布 ; Re=1200 时 , 产生 二 
次 回流 的 区 域 ， 温 度 发 生 了 波动 。 这 两 种 状态 的 共 
同 点 是 ， 温 度 在 台阶 高 度 以 上 通道 无 明显 变化 。 青 
附着 点 处 ， 流 体 冲 刷 壁 面 ， 壁 面 温 度 边界 层 破坏 ， 
温度 梯度 最 大 。Re 数 越 大 ， 冷 、 热 流体 之 间 的 换 热 
作用 增强 ， 两 者 温差 沿 着 流动 方向 减 小 得 越 快 。 

图 10 表示 温度 沿 通 道 站 法 向 的 变化 情况 。 顺 着 流 
动 方向 选择 9 个 截面 ， 其 中 截面 WH=0 位 于 台阶 分 
离 点 ， 截面 WEB=2、5 位 于 主 回 流 区 内 ， 截 面 xWH=7、 
10 靠近 再 附着 点 ， 其 余 截 面 位 于 主 回流 区 下 游 。 从 


高 度 方 


< 


Re 数 j 


i 


靠近 分 离 点 处 ， 温 度 只 在 台阶 高 度 以 下 
高 度 以 上 无 变化 。 主 回流 区 内 ， 温 


的 变化 高 度 逐 渐 扩 大 ， 台 阶 下 游 通道 上 部 
冷 流体 的 区 域 逐 渐 减 小 。 主 回流 区 下 游 ， 沿 通道 的 


， 流 体 的 温度 相对 于 进口 处 流体 温度 都 有 
所 提高 。 但 温度 在 靠近 底面 的 区 域 会 发 生 急 剧变 化 ， 
越 大 , 温度 的 变化 梯度 也 越 大 ， 随 着 高 度 增 加 ， 
昌 度 变化 越 来 越 平 组。 顺 着 流动 方向 ， 上 部 冷 流体 


和 下 部 热流 体 之 间 的 温差 逐渐 减 小 。 


3 结 


论 


采用 非 定常 方式 数值 模拟 了 二 维 后 向 台阶 内 部 


流 场 ， 


对 


台阶 底面 


的 传 热 特性 规律 和 通道 内 温度 分 


布 规律 进行 了 研究 ， 得 到 以 下 结论 : 

1) 在 层 流 状态 下 ， 台 阶 下 游 底面 的 瞬时 Nu 数 
和 时 均 Nun 沿 着 流动 方向 的 分 布 完 全 重合 。 

2) 在 过 渡 流 状态 下 ， 人 台阶 下 游 底 面 的 瞬时 Nu 


数 始终 围绕 时 均 Nunm 数 的 分 布 发 生 波动 ， 沿 着 流动 


方向 出 现 多 个 峰值 。Re 数 为 700 时 ， 主 


Ee 


流 区 未 出 


现 二 次 回流 ,瞬时 Nu 数 在 主 回 流 区 下 游 才 发 生 波动 ; 
Re 数 为 1200 时 ， 主 回流 区 出 现 二 次 回流 ， 了 瞬时 Nu 
数 在 相应 位 置 就 发 生 了 剧烈 波动 。 


化 呈 明 显 的 阶梯 式 。 在 同一 截面 ， 温 度 在 靠近 底面 


3) 后 向 台阶 通道 内 温度 在 各 截面 上 治 法 向 的 变 


处 发 生 急 剧变 化 ，Re 数 越 大 , 温度 变化 梯度 也 越 大 ， 
随 着 高 度 增加 ， 温 度 变化 越 来 越 平缓 。 
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